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With the appropriate choice of parameters and sufficient cooling, charged par- 
, Qh | tides in a circular accelerator are believed to undergo a transition to a highly- 

ordered crystalline state[l]. The simplest possible crystalline configuration is a 
one-dimensional chain of particles. In this paper, we write down the quantized 
version of its dynamics. We show that in a low-density limit, the dynamics is 
that of a theory of interacting phonons. There is an infinite sequence of n-phonon 
interaction terms, of which we write down the first orders, which involve phonon 
scattering and decay processes. The quantum formulation developed here can 
serve as a first step towards a quantum-mechanical treatment of the system at 

""J/5 ' finite temperatures. 
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1 Higher-Order Dynamics of the Coulomb Chain 



We consider an ensemble of charged point particles forced into an one-dimensional setup by 
an external focusing field. In equilibrium, the particles will be equidistant longitudinally. 
We treat the limit of an infinite, but periodic, chain. The problem will be treated in the rest 
frame of an orbiting particle, curvature and retardation effects will be neglected. 
The kinetic, potential and Coulomb interaction Lagrangian are, respectively 
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where we have introduced local coordinates around each particle's equilibrium position. The 
sums run over all lattice sites. A is the lattice vector, we use a coordinate system where 
A = (0,0, A). The particle has mass m, and the external focusing strengths are given by 
and are assumed to be constant along the ring. We are using natural units with 



ext,xi ^ext,y 

h = c=i. 

We expand (1) in that is, we write 
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df^ diverges, but is irrelevant here; Cf' = 0, as the coordinates are expanded around 
their equilibrium. For the first interesting orders, we get 
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where we used the shorthand notation A 1 = {8^ im — 5^ n ). 
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Doing the summation over m, n, we get 
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where $^ = |e 2 A n (sgn(yU — z/)) n 1 |/i — z/| n \ and = and II denotes all permuta- 
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tions of a set of indices. 
The sums over <3>fco give 
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for odd n and vanish for even n. (C(3) ~ 1.202, £(5) ~ 1.037). 

As the interaction is translationally invariant, we proceed by Fourier transformation: 
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We write down the interaction Lagrangian in this basis. For convenience, we introduce 
vertex functions 
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Note that momentum conservation is only up to integer multiples of 27T. After some Fourier 
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gymnastics, we have 
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2 Quantization 

The quadratic terms of the total Lagrangian describe an ensemble of harmonic oscilla- 
tors with coordinate variables ^(k) , Q (k) = We introduce momenta variables 
7Tj(fc),7r*(fc) = 7Tj(— fc) obeying the usual commutation relations. 

Quantization is straightforwardly done by defining creation and annihilation operators 
ai(k),af(k) by 

y/2n(k)aj(k) = n(k)£j(k) + iTTj(-k) , (9) 

with oscillator frequencies Q(k) defined below. These oscillator eigenmodes describe phononic 
(particle displacement waves) excitations of our system. 

We write the full Lagrangian (8) in terms of the operators Oi(k), at{k). The momentum- 
independent terms in (8) are disposed of by absorbing them into the Fourier transform of 
the potential: $(£;) -> $(fc) - $(0). 

Inspecting (8), one notices that the terms can be interpreted diagrammatically: 

1. gives the one-particle propagator, i.e., it gives the dispersion relation Vt 2 {k) for 
the phonons (Fig. 1) 

(3) 

2. F{ describes a decay process: one incoming phonon decays into two outgoing ones 
(Fig. 2) 

3. describes a decay process: one incoming phonon decays into three outgoing ones 
(Fig. 3) 

4. describes a scattering: two incoming phonons exchange momentum (Fig. 4) 

Note that our diagrams are in terms of the spatial coordinates £, £*. If we want to draw 
the diagrams in terms of phononic eigenmodes, we have to use £, £* oc a + ± a and draw all 
8 possible three-point and 32 possible four-point diagrams: Each leg in any of the diagrams 
can be flipped over to make an outgoing particle an ingoing one while changing the sign of 
its momentum. 
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Figure 1: Free two-point function 



Also, we have to multiply each diagram by the polarization tensors, i. e. the totally sym- 
metric ii dependent terms in (8). With an obvious notation for transverse and longitudinal 
polarizations, these are given in Table 1; contributions with index configurations not given 
in the table vanish. 



Index Structure 


Weight 


(-L.-L) 


-1 


(11,11) 


+2 


(II, II, II) 


-6 


(IU.-L) 


-2 


(-L,-L,-L,-L) 


+9 


(-L,-L,-L',-L') 


+3 


(-L.-L, IUI) 


-12 


(11,11,11,11) 


+24 



Table 1: Weight factors of different polarizations 



Looking at the coefficient of the two-particle diagram, we can write down the oscillator 
frequencies due to the internal degrees of freedom: 

2 00 

The explicit form of this dispersion relation involves generalized Zeta functions and is 
not too enlightening. However, we can write down the energy of the n mode, which is easily 
seen to be the highest energy mode (see Fig. 6): 
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Figure 2: Decay diagram 




Figure 3: Decay diagram 




Figure 4: Scattering diagram 
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Figure 5: Momentum insertion by the external lattice 




One reads off that £l\ jk = — and fi| k — 2Q|, that is, in the absence of external forces 
the transverse motion is unstable. This has to be counteracted by an external focusing field 
with a field gradient greater than f^.. In real- world situations, this field will be position- 
dependent, i.e., our Lagrangian ceases to be diagonal in the Fourier basis. Instead, we have 
a convolution with the Fourier decomposition of the lattice focusing. Diagrammatically, this 
means that the two-point functions can get injected momentum from the magnetic lattice 
(Fig. 5), the K ext (0) contribution just being the average focusing strength. 

Also, tt(k) determines the validity of our quantization procedure. Instead of quantizing 
the fermionic particles, we have quantized their collective phononic excitations, which we 
obtained by expanding the classical Lagrangian around the classical equilibrium (cf. [2, 3]). 
Obviously, the particles have to be localized even in the quantum-mechanical domain for 
this procedure to be valid. 

We have seen the particles behave oscillator-like to lowest nontrivial order. Thus, we 
can estimate their wave functions' longitudinal extension; the ground state of an harmonic 
oscillator has 

° 2 = — (12) 
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as its extension. For a point-particle expansion to be valid, we have to require 
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or, with (11) 
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(which is the one-dimensional version of the r s ^> 1 condition known from Wigner crystal 
theory [4, 5]), so the quantization procedure is valid for low particle densities. As the condition 
is expressed in the rest frame (so A = 7A^ a b), the condition can easily be fulfilled in realistic 
setups. 
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